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It is shown that the charged massive Schwinger model supports a periodic vacuum structure for 
arbitrary charge density, similar to the common crystalhne layout known in solid state physics. The 
dynamical origin of the inhomogeneity is identified in the framework of the bozonized model and in 
terms of the original fermionic variables. 
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^ , An important feature of Quantum Field Theories is the formal separation of the " active" and " spectator" degrees 
of freedom. In fact, the conventional contributions to the perturbation series of a Green function, represented by 
Feynman graphs, involve more and more actively participating particles as the order of the expansion is increased 
and the ground state with its infinitely many real or virtual particles remains formally hidden. This scheme is really 
efficient only if the vacuum is " empty" which is usually realized by the adiabatic turning on and off the interactions 
, as the time evolves. But serious complications arise in this scheme when the vacuum is "non-empty", i.e. contains a 
^ significant number of particles. When the constituents of the vacuum form a rigid system then the vacuum is called 
solid and space-symmetries are violated. When the constituents are not localized then the vacuum can be considered 
as liquid. 

When the fluctuations are sufficiently small then the saddle point expansion can be used to turn this qualitative 
picture into a systematical description. The "non-empty" vacuum consists of a condensate in this scheme. What 
is the momentum of the condensed particles? In case of vanishing momentum the vacuum is homogeneous and the 
> dynamics of the excitations can be described in a straightforward manner. But it may happen that the vacuum 
■ is made up by particles of non-vanishing momentum in which case the saddle point is inhomogeneous and breaks 
-j-J , the space-time symmetries. Depending on the 'inertia' of the saddle point the zero modes arising from this formal 

• symmetry breaking are either integrated over (liquid) or kept at a fixed value (solid). 
^ I A few examples for liquids are the following. The short range correlations of the vacuum in Yang-Mills models 
i correspond to a liquid of localized tree-level saddle points, instantons [1]. The one-loop effective action indicates 
' the presence of other condensates [2] which must be inhomogeneous in order to preserve the external and internal 
symmetries. The mixed phase at first order phase transitions and the corresponding Maxwell-cut results from the 
soft modes which are generated by the inhomogeneous saddle points of the Kadanoff- Wilson blocking procedure 
in renormalizing the action [3]. The metallic lattice is the best known example of solids, and similar, periodically 
modulated ground state belongs to the Wigncr lattice of dilute electron gas [4] or the charge density wave phase in 
solids [5]. 

Our current understanding of such vacua is severely limited due to the strong interactions or correlations between 
the "active" and the "spectator" particles. This is easy to understand in the framework of the saddle point expansion. 
The soft zero modes of the inhomogeneous saddle points in a liquid are easy to excite and they usually lead to serious 
IR divergences in the semiclassical expansion. There may not be soft modes present in the ground state of a solid but 
momenta borrowed from the inhomogeneous condensate generate nonperturbative phenomena such as the opening 
of gaps. In addition to the strong couphng between the "active" and "spectator" degrees of freedom the dynamical 
origin of the condensation and the explicit construction of the ground state from first principles represent a so far 
unmatched challenge in both cases. We believe that the treatment of the soft modes is a more serious and difficult 
problem than those of the momentum nonconserving processes. The goal of the present work is to trace the origin 
of the periodically modulated vacuum in one of the simplest interactive theory. Quantum Electrodynamics in 1+1 
dimensions, the Schwinger model [6]. 

Which part of the (effective) action is responsible of the inhomogeneity of the vacuum? The inhomogeneity suggests 
the presence of strongly distance-dependent interactions in the system and it is natural to expect that such interac- 
tions arise from the higher derivative terms in the action rather than from the ultralocal potential energy. Such a 
relation between the presence of higher order derivative terms with sufficiently strong coefiicients in the action and 
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the inhomogcncity of the vacuum has aheady been confirmed in a number of eases [7] . What was left open by these 
works is the dynamical origin of the higher order derivative terms in the effective action which drive the condensation 
at nonvanishing momentum. The higher order derivative terms which are supposed to be responsible for an eventual 
inhomogeneity in QED are to be sought in the effective theory for the photons or for the density. The simplifications 
which occur when we constrain ourselves into 1+1 dimensions allow us to use simple but powerful analytical and 
numerical methods, such as bosonization and the variational approach, to explore these effective theories in a nonper- 
turbative manner. The main result of this work is that the vacuum of the massive Schwinger model in the presence 
of nonvanishing total electric charge is periodically modulated. It is reassuring that both the bosonization and the 
variational approach yield the same conclusion. The bosonization allows to identify a mechanism, more involved as 
anticipated in the works [7], as the driving force in forming the periodically modulated vacuum. The inhomogeneity 
arises from the competition of an unusual piece of the kinetic energy which contains the first power of the space 
derivative, the boundary conditions and the periodic part of the potential energy. The periodicity of the vacuum field 
configuration in space originates from the periodicity of a term in the potential energy in the field variable. 

The Schwinger model has already been extensively investigated. What will be important from the point of view 
of the present work is that the confinement of the electric charge has been established [8-11] and simple analytic 
considerations hint that the vacuum with nonvanishing background charge is inhomogeneous [11]. 

It is worth mentioning that a non-trivial vacuum structure appears in QCD4 as well at high fermion densities within 
the large Nc expansion [12]. This result motivated the search for periodic structure in two-dimensional field theoretical 
models where the investigation can be carried out without approximations. It was found that the Gross-Neveu and 
the 't Hooft models exhibit periodic baryon density [13], and the multiflavor Schwinger model and QCD2 also shows 
up periodic ground state [14]. 

The resemblance of the Schwinger model to a one-dimensional electron system may one also lead to the idea of 
the existence of a periodic ground state. Some compounds can have such atomic structure that they exhibit one- 
dimensional metallic properties and show periodic structure [15]. An analogue of the Wigner crystal appears in a 
one-dimensional spin system with short range, nearest neighbour interaction [16], and a one-dimensional electron gas 
with a long range U{r) ~ 1/r type Coulomb potential also exhibits periodic ground state [17]. Another indication 
of the inhomogeneous vacuum structure in one spatial dimension comes from the non-relativistic Peierls mechanism 
[18]. 

The massless Schwinger model is exactly solvable [19] and the explicit computation of the fermion determinant 
leads to an effective theory with massive photons and confinement [9]. When the vacuum polarization effects are 
neglected then the electric flux conservation induces a flux tube between an electron-positron pair in the absence of 
other charges. The resulting linear potential, the hard confinement mechanism, renders the positronium confined. 
Let us now allow the vacuum polarization to be present and try to separate a member of the positronium, that of 
the meson of the Schwinger model. The electric flux tube breaks up due to electron-positron pair creation when the 
energy of the stretched flux tube is sufficiently large and the members of the newly created pair bind to those of 
the original pair. This is the soft confinement mechanism and one ends up with more neutral mesons again in this 
manner. The linear potential between the electron-positron pair becomes saturated by virtual pair creations. The 
potential between a pair of static test charges can easily be obtained in the presence of vacuum polarization [9] and 
it shows that the total screening, the soft confinement mechanism, occurs for arbitrary value of the charge. 

The massive Schwinger model is not exactly solvable and the potential between a pair of test charges is periodic 
function of the charge with period length given by the elementary charge e and is saturated for integer multiples of 
e only [9]. The mass gap prevents the vacuum to screen out non- integer multiples of the elementary charge. Notice 
that the massless model is singular in the sense that arbitrarily small mass is enough to prevent the screening of 
non-integer charges at sufficiently large distances. At the end any charge is confined in the massive model as well but 
the integer or non-integer charges are confined by the soft or hard mechanisms, respectively. 

Excitations above a fermionic vacuum with well-defined particle number are always consisted of particle-hole pairs 
and arc therefore of bosonic nature. It is the special feature of the 1+1 dimensional world that the effective theory for 
these bosonic excitations is local. The local effective theory resulting from the bosonization of the massive Thirring 
model is the sine-Gordon model [20] . These bosonization rules are widely used for the investigation of the Schwinger 
model [10,11,21-23]. The massless case yields a free scalar theory and the massive theory leads to the massive 
sine-Gordon model [24]. The massive Schwinger model was also investigated by bosonization technique [25] and by 
functional methods [26] at non-vanishing chemical potential and temperature T. The existence of a periodic chiral 
condensate with the wavenumber of 2/i has been established for arbitrary temperature, too. 

Arguments were given in the framework of the tree- level solution of the bosonized model [11] that the massive 
Schwinger model exhibits a periodic ground state in the presence of a static, homogeneous background charge density. 
A more systematic investigation of the inhomogeneity of the vacuum of the massive Schwinger model in the presence 
of homogeneous external charge density, Pext) is presented in this paper. We attack the problem from two different 
directions. First by minimizing the tree-level expression of the energy functional in the bosonized form of the model 
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and after that by minimizing the energy with respect to the parameters of a static, periodic background electric 
field in the fcrmionic form of the theory and by retaining the quantum fluctuations of the photon field up to the 
two-loop order. The results obtained by both approaches are in qualitative agreement. Namely, the ground state 
exhibits periodically modulated charge density with decreasing amplitude for increasing pext- For large Pext numerical 
calculations failed to be conclusive regarding the true energy minimum. Analytic considerations were used in this 
density regime with the result that the vacuum remains periodic for arbitrarily large values of pext • 

The paper is organized as follows. Section II contains the study of the tree-level bosonized theory in the presence 
of the homogeneous external charge density poxt • The minimum of the energy functional of the model is found by 
numerical minimization of the classical vacuum energy. For large values of poxt when the result is more unstable 
with respect to numerical errors the stability of the periodic vacuum has been shown analytically, by expanding the 
tree-level energy in powers of the amplitude of charge density wave in the vacuum. 

The fermionic investigations are presented in Section III. Since integer charges are screened and do not leave behind 
long range photon field [11] the perturbation expansion in e is reliable by using the original fermionic and photonic 
degrees of freedom. We follow a variational strategy and minimize the energy of the vacuum as the function of 
the induced photon field in the vacuum. The fermionic degrees of freedom are integrated out in the presence of a 
static, sinusoidal electric field and the quantum fluctuations of the photon field are taken into account up to two-loop 
diagrams for the energy. The external charge density pext is introduced indirectly via a fermionic chemical potential 
/i. The energy of the vacuum is finally minimized with respect to the amplitude and the wavelength of the static 
periodic background electric field. The numerical minimization procedure finds the periodic ground state energetically 
favorable as compared to the homogeneous one below certain value of Poxt- The problem of the high density regime 
is of the same origin as in the bosonized study, namely the smallness of the amplitude of the induced photon field in 
the vacuum. An analytic calculation in the framework of the perturbation expansion in the amplitude of the induced 
photon field predicts a periodic ground state even in this density region. 

Finally, the conclusion is drawn up in Sect. IV. Appendices A and B present briefly the numerical solutions of 
the Dirac equation in periodic background potential and the band structure of the fermionic spectrum, respectively. 
Explicit expressions for the Feynman diagrams for the energy and charge densities up to the two-loop order are given 
in Appendix C, and the details of the numerical search for the energy minimum are given. 

II. MINIMIZATION OF THE ENERGY IN THE BOSONIZED MODEL 

This section contains the tree-level determination of the vacuum structure of the bosonized model. 



A. Hamiltonian 

The Lagrangean of the massive Schwinger model is given as 

C = -\f^.F'^-' + i>r{d^ - ■mA^)i, - m^V, (1) 

where F^i, = d^^Ai, — d„A^, m and e are the bare rest mass of the electron and the bare coupling constant, respectively. 
The bosonization rules are [11]: 

: tptjj : — cmM cos{2^/7r<j)), : xp^^tp : ^ — cmM sin(2-\/7r0), 
=■ V^l/.V' : ^ -^e^^d'^cj), : V^i^ : ^ \N^{d^<f>)\ (2) 

where denotes normal ordering with respect to the fermion mass m, c = exp (7) /27r with the Euler constant 7, 
and M = c/i/tt the 'meson' mass. It is believed that the presence of a non- vanishing background charge density does 
not affect these transformation rules [25]. The Hamiltonian of the system in Coulomb gauge is given by 

W = / i)x(yiidx + m)i>x - ^ ( jQ,x\x-y\jo,y, (3) 
Jx 

with = dx^ J^j^ dx^ . According to the bosonization rules this Hamiltonian is equivalent to those of the massive 
sine-Gordon model, 
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n[ii,(f>\ = N„ 



:{di(l)x) 



-M'^(f>l - cmM cos{2s/n^x 



(4) 



where 11^: denotes the momentum variable canonically conjugated to 

Our purpose is to determine the vacuum of the massive Schwinger model in the presence of an external static 
particle density Pext x which is added to the density jo,x in Eq. (3), 
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0,x + Pext x)\x - y\{jo,y + Pext y)- 



The external charge is represented by the external field (^ext x in the bosonized Hamiltonian (4) as 



J X 



hil + ^{di(t>xf + ^M^i^x + '/'ext xf - cmM cos{2^(j)x) 



where 



(5) 



(6) 



Pext 



di(, 



(7) 



The external particle density is assumed to be static and constant in the interval = z £ Y^L, L] and vanishing 
elsewhere, therefore we write </)ext x = bz when |z| < L and </)ext a = elsewhere for any =t. It is advantegous to 
introduce the field variable 



which allows us to write the total particle density as 

Px = 



(8) 



(9) 



The tree-level vacuum can be constructed by minimizing the Hamiltonian (4) as the functional of the static field 

configuration (j)x with n^; = 0. The minimum is reached at c^gr x = (0|</>,i|0) and the value of the Hamiltonian at this 
field configuration, E(b) = ?iext[0, (/igr], can be identified by the tree-level vacuum energy. Lattice regularization of 
the Hamiltonian for static field, 
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J z 
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TiiMzf + 7r^4>z + bzf - COs(2v^0^) 



yields 



N 2 ^ ^ 

aHLiO, 4>] = - Y.{4>n+1 - + ^ Y.^4>n + hZnf - ^ C0s(2^AF.^„), 

n=0 ra=0 ^ n=0 



(10) 



(11) 



where a stands for the lattice spacing, cl = ea, tul = ma, Zn = zq + na, a = 2L/ {N + 1), zq = —L, Zn+i = L and 
•^n = 4>zn ■ The boundary conditions 



4>u = <Pn+i = 



(12) 



have been used in order to restrict the computation into the sector with vanishing induced charge. 

In order to understand the origin of the periodic structure of the vacuum we rewrite the static Hamiltonian (11) in 
terms of the shifted variable 6 as 



N 



\ X^(0n+1 - - V^Pext i)^ + ^ " ^^^^^ 

n=0 n=0 ^ " 71=0 



N 



* r» — 



cos(2V7r(^„ - hz)), 



(13) 



where pext l = apext = ahj \fv: denotes the amount of external particles distributed between two consecutive lattice 
sites. This expression reveals a competition in forming the vaciium, taking place between the kinetic and the potential 
energies, the first and the remaining terms on the right hand side of Eq. (13). Let us first ignore for simplicity the 
quadratic mass term and the shift —hz in the argument of the cosine function on the right hand side, a simplification 
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which yields the Haniiltonian of the Frenkel-Kontorova model [27]. The tree-level vacuum of this model produces 
infinitely many commensurate-incommensurate transitions and displays a rather involved phase structure with the 
devil-staircase feature due to the competition between two dimensionless parameters [28]. In fact, the kinetic energy 
prefers 

= const. + n\/7rpext L (14) 

and the potential energy is minimal for <pn = J^/tt (with n and j integers) and the vacuum is trivial, i.e. a linear 
function of the coordinate, for integer pext L only. The vacuum of the complete Hamiltonian (13) is always the result 
of a compromise between the kinetic energy with the preference expressed by Eq. (14) and the potential energy which 
prefers (6 = and ^„ = \/7r(j — npext l)- The competition between the kinetic and the potential energy is never trivial 
due to the quadratic mass term but is at least simplified when pext l is integer. For sufficiently small lattice spacing 
Pext L < 1 each energy expression on the right hand side of Eq. (13) enters in the competition for the vacuum. 

Such an involved vacuum structure is characteristic of the tree-level solution in lattice regularization only. The 
quantum fluctuations should smear most of the commensurate-incommensurate transitions out. A similarly smeared 
behaviour is what one finds when the cutoff is ignored in the tree-level sector, i.e. when the minimum energy 
configuration is searched in the naive, classical continuum limit of Eq. (13), 
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t[0,^]= / 

J Z 



Adi(l)z - V^r/Oext)^ + TT'i^l 1= cos(2V7r(02 - ^/i^Pe^tz)) 



(15) 



subject of the boundary conditions ^-l = ^, 4>l = Zy^-Z^Pext- In the Frenkel-Kontorova limit when the mass term 
and the ^-dependence in the argument of the cosine function are ignored then (j)^ develops oscillatory structure in the 
vacuum which changes smoothly with the parameters of the model. In fact, the kinetic energy prefers to distribute the 
total change 4>l — 4>-l = 2-\/7ripext in a linear manner but the potential energy introduces a periodic modulation. The 
period length can be determined by noting that (j) should change by within a period. Such a simple argument gives 
the period length £o = l/Pext for small me. The fermions correspond to kinks of the sine-Gordon model according to 
the bosonization therefore it is not surprising to find that there is just one particle per period in such a vacuum state. 
The z-dependent shift in the cosine function takes out the driving linear term from but tends to generate periodic 
oscillations with the same period length. 

Notice that the source of the inhomogeneity of the vacimm is an unusual. 0{di), gradient term in the kinetic energy. 
This contribution to the energy, together with the boundary conditions and the periodic potential energy form the 
periodic modulation in the vacuum. Furthermore the space inversion symmetry is broken explicitly by the 0{di) term 
and the boundary conditions. It is interesting to compare this situation with those encountered in earlier studies [7] 
where the dispersion relation of the form 



(p) = C6/ + C4P^ + ^+Co (16) 



was used with Ce, — C4 > with non-periodic potential energy and the tree-level vacuum is expected to be periodic 

if there is a region in the momentum space where e(p) < 0. The space-time inhomogeneitics are therefore generated 
by the competition of terms with different orders of the gradient only and the space inversion symmetry is broken 
spontaneously. 



B. Numerical results 



The energy minimum was searched by the conjugate gradient method which started from a number of initial 
conditions for (pzi and the field configurations corresponding to the lowest energy only have been singled retained. 
The charge densities were then calculated according to Eq. (9). We used L = 167r, A'' = 800, m = 0.5; 2; 5 and 
b G [0.3; 7] in the numerical studies. 

The minimum of the expression (10) was found at = for vanishing external charge density, i.e. for 6 = 0. The 
increase of b gave two distinct regions, separated by a size-dependent point b = 6^. 

b <br. : For b close to zero one expects that <pz is small and the sinusoidal potential in the equation of motion, 

di(t)z = —i^z + bz) + 2cmesm{2y/^(j);,), (17) 
can well be approximated by the first term of its Taylor series 
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e2 



di(j)z ~ —{<Pz + bz) + 4cme\/7r(?f)2. (18) 

TT 

Such a linearized equation of motion together with the boundary conditions (12) yields 



sinh(Kz) I ^ — ^ 

4>z=0sL . - bsZ, K=\ l-4cmeV7r, bs = b^^. (19) 

Sinh(KL) V TT K^TT 

This solution, shown in Fig. 1, contains three spatial regions in the interval [—L, L] and in the longest, central region 
(f)z is linearly decreasing function with the slope —6s. The analytic results for the slope are in very good agreement 
with those obtained numerically. The linear decrease of in the central region describes the partial screening of the 
external charge density. In the two other regions, close to the boundaries at —L and L, \(j)z \ approaches abruptly zero. 
For every choice of i a critical 6i, value was found where the linear approximation fails to work and the higher-order 
terms of the sine function arc needed in the equation of motion (17). It was found that the slope reaches b at this 
point. The L-dependence is 6l » L~^-'^^ according to Fig. 2, therefore 6l ^ and this type of solutions disappears 
in the thermodynamic limit. 

br, < b : The increase of bg to b indicates the complete screening of the external charge density in the central region. 
Furthermore the numerical solution, depicted in Fig. 3, reveals an additional periodic structure in (pz, 4>z is a periodic 
function of wavelength £, (j)z = 4'z+t- The wavelength (. and the amplitude A of were defined numerically as the 
distance of the neighbouring zeros of (f)z and the arithmetic average of the magnitude \<f>z\ at the extrema of the 
periodic component, respectively. Both i and A decrease with increasing 6 in this region as shown in Figs. 4 and 5. 
This feature opens the possibility of applying the perturbation expansion in the amplitude A of the induced periodic 
field in the vacuum in the limit of asymptotically large charge densities, p — > oo. Based on Fig. 3 the periodic part 
of the scalar field is approximated by 

^, = ^sinf^^V (20) 



By inserting this expression into Eq. (10) one finds the energy density 

where Ji(x) is the Bessel function of the first kind. Fig. 4 shows that the amplitude A decreases with increasing 
charge density. Therefore it is sufficient to consider the expression on the right hand side of Eq. (21) only up to the 
order 0{A'^) for large p ^ oo. Due to the relation 

J,{x)^^ + 0{x^) (22) 

valid for small x the energy density takes the form 

S{A,b) = ^ + ^ f 47r62 + —)- cmeA (23) 
2 4 V TT J 



having a non-trivial minimum at 



2cme „ 
47r62 -I- £- 



where the Casimir energy is negative, 

E{A,b)-£{A = 0,b) = --^—^<0. (25) 

TT 

Thus one concludes that the ground state of the massive Schwinger model is periodic for large external charge 
densities. It is shown in Fig. 4 that the analytic result of Eq. (24) is in good agreement with the numerical one for 
the charge-dependence of the amplitude A. 

The relation £ = 1/pext displayed by Fig. 5 reflects the fact that charges which are integer multiples of e are 
completely screened. In fact, as argued in Ref. [11], the introduction of the charges ±e at the boundaries corresponds 
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to the shift z ^ z + 5z with \Sz\ = y^/b. This is a symmetry of the vacuum if has the length of period £ = \Ez\lv 
where v is integer. According to the numerical results i/ = 1 at the energy minimum. Similar periodic structure is 
found in Wigner crystals of itinerant electrons, in certain spin systems [16] and in the charge density wave states. 
The periodicity usually gives way to homogeneity when the external charge density is increased because the overlap 
integrals between the neighbouring lattice sites increase. This is not what happens in the massive Schwinger model, 
where the simple, leading order perturbation expansion given above shows that the ground state keeps its periodicity 
for arbitrarily large charge densities. 

Our conclusion is that in the tree-level approximation of the bosonized theory the massive Schwinger model has 
a single periodic phase in the thermodynamic limit and the homogeneous external charge density is neutralized in 
average by a periodic, induced charge density. Integer charges are completely screened as argued in [11]. 



Let us consider now the massive Schwinger model in terms of the original fermionic degrees of freedom and subject to 
periodic boundary conditions at the endpoints of a finite spatial interval. The finite charge density is now introduced 
by the chemical potential /z. The system of electrons is easier to polarize than the 'empty' vacuum and accordingly 
there is no gap in the free electron excitation spectrum for /i > m. The photon polarization tensor is non-vanishing at 
the Fermi level therefore the Debye screening renders the photon propagator short ranged and the Coulomb potential 
vanishing for large separation [29]. Our computation performed in this formalism supports the results obtained in the 
bosonized theory, namely that even an arbitrarily weak interaction among the electrons is sufficient to form a periodic 
ground state. The dynamical origin of the modulated ground state is the opening of a gap around the Fermi level. 

Wc are confronted by two complications in describing the vacuum. First, the confinement of charge renders the 
fermionic excitation spectrum non-physical and ill-defined. As discussed above in the framework of the bosonized 
theory integer multiples of the elementary charge are screened by vacuum polarization at finite charge densities and 
their Green function is short ranged. Since only integer charges can be created in the fermionic theory wc expect no 
problems with perturbation expansion at finite density. The second problem, the possibility of dynamical generation 
of coherent photons, i.e. a background field in the vacuum is more difficult and has to be handled in a self-consistent 
manner. For this end we introduce an external photon field. 



with a > chosen to be a single plane wave for the sake of simplicity. Since there is only one non- vanishing component 
of the field strength tensor Ff^^, = d^Ai, — d,^A^, such a background field represents a generic sinusoidal external field in 
1 -|- 1 dimensions. The energy density will be computed in the order 0{e^) in the vacuum for a given /i and minimized 
with respect to the external field, the variational parameters a and Q. 

The numerical minimization of the vacuum energy density with respect to the background field shows that the 
system manages to lower the vacuum energy below the 'empty', perturbative value by opening a gap and generating 
a photon condensate (26) in the vacimm for small densities. For large densities the perturbative treatment of the 
dependence of the vacuum energy on the field (26) is reliable and yields similar results. Our analysis does not cover 
the intermediate density regime where the density is large enough to make the numerical minimization of the two-loop 
energy expression unreliable but small for the application of the perturbation expansion in a. 



The background field is introduced by the collective coordinate method into the generating functional for the Green 
functions. The vacuum-to-vacuum amplitude of the model is expressed by the path integral 



III. VARIATIONAL MINIMIZATION OF THE ENERGY FOR QED 



(26) 



A. Background field as collective coordinate 




(27) 



where the action for the photon field A'* in Feynman gauge, 




(28) 



is expressed in terms of the inverse of the free photon propagator 



7 



(i?-i)S^ = 5''^n.5.,„ (29) 

and the Dirac action, 

Sj,[A,ij,^] = ij-G-\A)-ij, (30) 
is given by means of the inverse fermion propagator 

G-\A) = ij^id^ - ieA^) -m = ^""ido - Hn{A) (31) 

with the Dirac Hamiltonian 

HD{A)=j°{-ij^di+m-erA^^). (32) 

We use the notation = dx^ dx^ ^ f ' 9 = fx9x and shall consider the limit LT oo below. 

The vacuum of the model will be constructed by means of a variational method. We introduce an external back- 
ground field AJ^ and separate the qiiantum fluctuations a'^., A^^. = A'^^ + a'^. The dcpcmdcnce on the external field is 
retained by the method of collective coordinates which implies the insertion of the identity 

1 = j daS{C[A, a] + a) (33) 

into the path integral, 

Z = J daZ„, j ^[V']2^[V']2?[a]'5(C[^, a] + (j)ei'SEM[A+a]+iSD[A+a,v;,V] (34) 



with 



C[A,a] = h^F-F)-F=^a-D-'^-A. (35) 



The fluctuations of the collective coordinate a are suppressed in the thermodynamic limit because the background 

field is extended and the cr-intcgration can be performed by expanding In Z„ around its maximum. The contribution of 
the collective coordinate to the vacuum energy density will be negligible in the thermodynamic limit and the collective 
coordinate can be frozen at the maximum as far as the energy density in the vacuum is concerned. 

One usually employs the effective action formalism in similar problems. There the external source, coupled linearly 
to the fluctuating fleld is supposed to stabilize the vacuum with the desired condensate. The minimization of the 
effective action guarantees that the external source plays no role in the true vacuum. The complication which renders 
this method rather involved beyond the leading order of the loop expansion is the Legendre transformation. The 
procedure outlined above leads to simpler expressions in the two-loop order. Both methods are useful in the case of 
stable ground state only. Large amplitude fluctuations appear in the mixed phase which make the computation of 
the convex effective action and the taking into account the fluctuations of the collective coordinate difficult. 

It will be useful to introduce the generating functional 

Za[j,C,C]= J dX J p[^]X)[V']X>[a]e'^™[^+"l+^^°[^+"''^''''l+^^(<='[^'"l+'^)+^^'-^ (36) 

where the constraint is represented as a Fourier integral over A. The generating functional can be written in the 
perturbation expansion as 

00 / \ " 

^>i*.<.<i= Ej('»i4^:44j ^-b-.c,ci, (37) 



where 



Zo a[j,C,C] = exp 



Tr InG-(A) - iC • G(A) • C - - (21 + 4i^) 



-^TrlnD-i - 1 In (-^AA) + ^Q^A-A - ^ ■ D' ■ j 



(38) 
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with 

The photon propagator (39) tends to the free photon propagator in the thermodynamic limit when the fluctua- 
tions parallel to the background field vanish and we continue using the original photon propagator D. Finally, the 
expectation value of an operator 0[A] is determined as 



0[A] = ^-^T^ JdxJ V[mmc^]0[A]e' 



SEM[A+a]+iSDlA+a,iP,i)]+i\{ClA,a]+a)+ija+iCtl)+itpC 



, (40) 

j=C=C=<r=0 



where the field variables in the operator 0[A] are replaced by functional derivatives with respect to the corresponding 
external sources. 



B. Energy and chcirge of the vacuum 

It may happen that the system prefers energetically a periodic ground state rather than the normal, homogeneous 
one. Then it should adjust itself by building up a static, periodic electric field and the corresponding band structure 
with the Fermi-level placed in a forbidden band. In order to decide whether such a readjustment of the vacuum takes 
place one has to compare the energy densities of the homogeneous and modulated vacua. The fermion spectrum is 
needed for the determination of the energy and charge densities. It has been calculated along with the corresponding 
eigenspinors by solving the Dirac equation numerically (see App. A and also App. B for a detailed discussion of the 
fermion spectrum). The fermion Green functions were constructed according to their Lehmann expansion. The poles 
on the complex energy plane were shifted according to the rules established in [30] to take into account the chemical 
potential ^. The energy density £[A] of the system is given by the expectation value of the 00th component of the 
energy-momentum tensor. At the two-loop order this expectation value can be represented diagrammatically as 

£[A] = ^a^Q^ + (^^^^ ~ ^ ^ i ^^^^^ ^ t^^^y..^^^ _ ^I^Sl^ ^ (41) 

the details of the calculation are given in App. C 1. There is no need of mass and charge renormalizations in QEDi+i, 
but UV and IR divergences appear in the diagrams, which should be properly handled, see App. C 2. We introduce 

the Casimir energy which is usually the energy difference of the states with and without a classical object and was 
investigated thoroughly [31] by the collective coordinate method. In our case it is the background field which plays 
the role of the classical object and the Casimir energy is 

£c(a,Q,M) =fpcr(a,O,/i)-f„(0,0,Ai). (42) 

In order to understand the structure of the vacuum, we need another important observable, the average charge density 
p. Its two-loop order expectation value is 





p[A] = ) + "{J) — {J) - ^wwvU , (43) 



for more details see App. C3. Charge neutrality for ^ = implies the renormalization condition Prenf^ = 0] = 
which can be satisfied by the subtraction 

p,en[A]=p[A]-p[0]. (44) 



C. Numerical results 



The energy densities given by Eqs. (41) and (42) have been calculated numerically for both the periodic and the 
homogeneous phases using the explicit formulae of App. C4 in units of e = 1 and the details of the numerical 
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procedure are discussed in App. C5. The negative value of the Casimir energy, found numerically, indicates that the 
periodic state is energetically favoured. The one- and two-loop contributions to the energy are shown in Fig. 6. The 
one-loop contribution, i.e. that of the first diagram on the r.h.s. of Eq. (41) is more important than the two-loop 
terms in the parameter range studied. The two-loop correction which is dominated by the exchange diagrams, the 
third and the fifth ones on the r.h.s of Eq. (41), tends to destabilize the periodic state for e/x < m and to stabilize it 
for e/i > m. The periodic phase is stable mainly due to the gain arising from the sinking of upper bands to negative 
energies and this gain is taken into account completely in our computation since the one-particle energy levels have 
been calculated non-perturbatively. The photon exchange appearing at the two-loop order becomes important in the 
region e/i « m. The numerical results for efi/m 1.12 agree for the periodic and the homogeneous ground states 
within the numerical accuracy. 

The computation of the charge density yields the equation of state, the relation between the energy and the charge 
density. For each /x and m wc looked for the values of a{^) and Q(/i) at which the Casimir-encrgy density assumes its 
minimum. In this way we obtained fc(a(M)) Q(a*)) m) = ^c(a*), and Pren{fi) for each fixed value of m. The equation 
of state then can be constructed by tracing fper(/u) and £n{l^) as the function of Pren(M) for the modulated and the 
homogeneous phases, respectively. A typical p-dependence is shown in Fig. 7 for m = 2, similar curves are found for 
the other values m. 

The two-loop results shown in Fig. 7a indicate in a reliable manner that the energy minimum for the periodic 

groimd state, fpcr, is smaller than that for the homogeneous groimd state, for p G [0,0.33]. The amplitude of 
the periodic background field tends to zero with increasing p as one can see in Fig. 7b. Our numerical results for 

p ~ 0.33 show rapid oscillations in Figs. 7b and c reflecting a numerical problem which arises due to the almost 

degeneracy found for the periodic states with Q = kp and 2kp. It was found that the wavenumber Q of the periodic 
phase is directly related to the Fermi momentum kp via the relation fQ — 2kp, with the filling factor /, defined as 
the number of the entirely filled positive-energy bands plus the fraction of the partially filled band, see Fig. 8. The 
ratio Q/kp displays a discrete behaviour, Q/kp » 2 for e/x < m and Q/kp ~ 1 for e/x > m. The discrete nature of 
Q/kp reveals that the vacuum always readjusts itself until a forbidden zone is opened at the Fermi-level and the filling 
factor becomes integer. This observation is in agreement with the nesting relation fQ = 2kp. The product fQ shown 
in Fig. 7c increases in the average with increasing charge density, implying the same behaviour of the wavelength as 
was found in the bosonized theory, c.f. Fig. 5. 

For densities larger than the 'critical density' p > 0.33 wc found no energy difference between the periodic and 
the normal ground state within our numerical accuracy. The reliable numerical determination of the parameters a 
and Q became impossible, so that those are not shown in Figs. 7b, c for large densities. The result obtained in the 
framework of the bosonized model, namely the small but non-vanishing amplitude of the modulation of the vacuum 
for large densities as shown in Fig. 4 makes one cautious that above the 'critical density' the numerical minimization 
fails to find the proper minimum. But fortunately the smallness of the amplitude of the induced field at large densities 
enables one to treat the periodic potential as a perturbation. The external potential (26) should open a gap around 
the Fermi level at kp = Q/2 according to the leading order of the degenerate perturbation expansion [32] and one 
finds the energy spectrum 



.(0) 



(45) 



close to the quasi- momentum k ^ kp, where e^"^ refers to the unperturbed spectrum and the signs -|- and — stand for 
|A;| > and |fc| < ^Q, respectively. Due to the first order perturbation treatment we can now ignore the two-loop 
order diagrams in Eq. (41). By means of the energy eigenvalues in Eq. (45) one finds 





dfc' 



Q/2 



271- hQ-t'' 



Q/2 



Q/2 



dk' 



2ti hQ-k" 



(46) 



for the Casimir-energy density in the leading order. Since p = /^^ dk/2'K = kp/iv and Q/2 = kp = wp we have 

p ~ Q/2'K. For large values of the chemical potential p wc expect large values of kp = Q/2 ^ m therefore we can set 
m = 0, i.e. e^''^(/c) = \k\, p, = kp = Q/2 and one can perform the integral in Eq. (46) explicitly. 



£c{a, Q, Q/2) = ^ + Q — z 



e^a^ . ^ Q 
— — arsmh - — . 
2tt 2ea 



(47) 



For a given Q, i.e. p the Casimir-energy density has two extrema: a maximum at a = {d^£c{0,Q,Q/2) — oo) 
and a minimum at 
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(c?^fc(«min, Q, (5/2) > 0). The latter provides the Casimir-cncrgy density of the ground state, 



^c(amin,Q,(3/2) = ^ 1 




2 " 

= — 1 - coth 
2 




< 0. 



(49) 



One finds flmin ~ 0.004 and 5c(amin, 27rp, 7rp) w — 10~^ at the 'critical density' p = 0.33 therefore it is not possible 
to confirm the periodicity of the ground state within our numerical accuracy. Furthermore we have seen for lower 
densities that the jump from one band to two bands in the Dirac sea is a two-loop effect caused by photon exchange. 
In general, one needs higher-loop corrections in order to let more than two bands sinking into the Dirac sea. This is 
due to the observation that the N-th order perturbation expansion in the monochromatic external field predicts the 
opening of N gaps in the fermion spectrum. Nevertheless this simple computation indicates that the vacuum of the 

Q 3 2/2 

massive Schwinger model keeps its periodicity, although with decreasing amplitude a ^ pe^"''^ p and wavelength 
I-kIQ = ^/ P for increasing charge density p, in a manner similar to the bosonized model (sec Fig. 5). The one-loop 
level vacuum for the assumed simple sinusoidal background potential involves a single band sunk into the Dirac sea. 

For the low-density periodic phase C/Lt « m where the amplitude a is non-perturbative the photon exchange is sig- 
nificant. As p increases the perturbative region of a is reached, although we could not decide, whether the numerically 
found 'critical density' does fall into the perturbative region. Our one-loop perturbative result cannot clarify whether 
higher-order loop corrections can lead to sinking more and more bands into the Dirac sea as the density p increases. 



The ground state for the massive Schwinger model has been investigated in the presence of homogeneous external 
charge density. The energy density of the ground state has been determined numerically in the bosonized version of 
the model as well as in terms of the original degrees of freedom of QED by variational methods. The scalar field 
configuration and the charge density of the ground state for the bosonized model have been obtained by minimizing the 
tree-level energy in the presence of static, homogeneous external charge density. In the fermionic theory a variational 
method has been constructed by minimizing for the amplitude and the wavelength of a sinusoidal photon condensate 
in the vacuum. The finite charge density was realized by the introduction of the chemical potential. The applicability 
of the loop-expansion for the computation of the vacuum energy at finite density is justified by the bosonized theory 
which indicates that the confining Coulomb-force among integer charges is vanishing. 

The computation in the bosonized theory shows that the system exhibits a periodic ground state for arbitrary charge 
density. The fermionic computation gives the same result. Numerical computations reliable up to a certain charge 
density as well as simple analytic calculations valid for asymptotically large charge densities support the periodicity 
of the ground state. The general trends in the charge-density dependence of the amplitude and that of the wavelength 
of the periodic structure are in agreement for both versions of the model. 

The investigations in terms of the bosonic and fermionic degrees of freedom complement each other. On the one 
hand, our results for the bosonized model showed that the background charge density is in average neutralized in the 
ground state. Furthermore, the charge density wave ground state and the complete screening of the integer charges 
appear due to the interplay of the kinetic energy, the periodic potential energy and that of the boundary condition 
for the boson field. The fermionic description, on the other hand, gives more insight into the structure of the vacuum, 
namely that the modulation of the charge density arises as the result of the opening of a gap in the fermion spectrum. 
The conclusion of the straightforward perturbation expansion for large charge densities is that the periodic ground 
state is always favoured by the system energetically as compared to the homogeneous one, even if only a single band 
sinks into the Dirac sea. In the case of a single mode periodic potential bands sink below the Fermi-level in the 
N-th. order of the perturbation expansion. In particular, in our numerical computation where a single mode was 
allowed for the induced photon field in the vacuum and first (second) order perturbation expression for the vacuum 
energy was minimized one (two) bands are found below the Fermi-level for large densities. The Peierls mechanism is 
present in the relativistic vacuum for N > 2. It is however not clear how to identify this mechanism in terms of the 
bosonic excitations. 

Our analytic considerations showing the existence of the periodic ground state for arbitrarily large charge densities 
do not include the corresponding higher-order loop corrections neither for the bosonized model, nor for the fermionic 
one. It still remains an open question whether a complete resummation of these higher-order corrections would alter 
the qualitative result obtained here, namely that the massive Schwinger model has only a single phase, the periodic 
one for any values of the average charge density. 



IV. SUMMARY 
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APPENDIX A: DIRAC EQUATION WITH SINUSOIDAL POTENTIAL 

The Lehmann representation of the non-interacting electron propagator requires the knowledge of the eigenfunctions 
of the Dirac Hamiltonian (32) for the external field (26). We use the real Dirac matrices 

The eigenspinors f^"{x) and g^^{x) of the Dirac Hamiltonian HoiA) belonging to the positive and negative energy 
eigenvalues efc^^>0 and — e^~^<0, respectively satisfy the equations 

Hn{A) f^\x) = eit^ f^x), Hn{A) g'^{x) = -e^;) g'^^x), (A2) 

where the quasi-momentum k G [—(3/2, Q/2) takes values in the first Brillouin-zone. The non-negative integer s > 
labels the bands in increasing order in the energy. As in the non-relativistic case, the solutions of Eq. (A2) are 
Bloch- waves, 

n=— OO n=— oo 

with kn=k-\-nQ. In order to find the numerical solution, one rewrites Eqs. (A2) in matrix forms for the components 
of the Bloch-waves, e.g. the first one of Eqs. (A2) reads as 

oo 

^ [{eks + e/x - fc„7%' - m7°) <^ + y + e-'^'"^^"-''"^^^ = 0. (A4) 

n=— oo 

The solution is found by making up a matrix from the coefficients appearing next to the Dirac spinors w^*. The non- 
relativistic treatment results in a matrix with tridiagonal structure [33], [34]. The structure of the matrix remains 
unchanged in the relativistic case except of the replacement of the matrix elements with 2x2 matrices. The problem 
is then reduced to solve a system of a coupled set of homogeneous linear equations. 



APPENDIX B: BAND STRUCTURE 



In order to understand the effects of the photon exchanges, the salient features of non-interacting electrons in static 
periodic background field are briefly summarized in this section. The eigenspinors and the energy eigenvalues of 
the Dirac Hamiltonian (32) in the static, periodic external field have been determined numerically (App. A). The 

positive and negative single-particle energies are denoted by e^"^^ and — e^~\ respectively, as the functions of the 
quasi-momentum k and the band index s. Their dependences on the parameters a, Q are not indicated explicitly. 
Since the chemical potential fj, results in a constant shift of the whole fermion spectrum, it is sufficient to understand 
the spectrum for /x = 0. Due to the periodic potential, a band structure with alternating allowed and forbidden bands 
is formed [35]. The typical band structure is plotted in Fig. 9a and 9b as the function of 1/Q for undercritical ea < m 
and overcritical ea > m static periodic external electric fields, respectively. The shaded regions in Figs. 9 correspond 
to the allowed bands. The non-relativistic analogue of these figures can be found in [33] where the one-dimensional 
electron system was considered in the presence of the static external electric potential Aq = a(l — cosQx) (with our 
notations) that is just the same potential we have but one of its minima is shifted to a; = and the potential is chosen 
zero in this minimum. The following qualitative features of the fermion spectrum are wothwhile mentioning. 
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1. Undercritical vacuum ea<m: The mass gap around zero energy separates the infinite towers of bands above and 
below this gap. The smaller is 1/Q, the more the allowed bands widen out and start to overlap. This is just the 
qualitative behaviour obtained in the non-relativistic description [33]. The energies of the states in the upper, 
e > (lower, e < 0) tower decrease (increase) with increasing 1/Q. 

2. Overcritical vacuum ea > to: The qualitative features described above remain the same but the upper and the 
lower bands overlap for large enough 1/Q and upper bands sink into the Dirac sea while lower bands emerge. 
The band crossing is well understood [29] but this case which involvs the creation of electron-positron pairs 
turned out to be not relevant for our purpose because the periodic ground state was found undercritical. 

The asymptotics of the spectrum for l/Q — > and l/Q — > oo is helpful to understand the l/Q-dependence of the 

band structure. 

1. 1/Q ^ 0: The energy levels decouple from the periodic structure of the potential and the free fermion spectrum 
reappears with the single gap for —m < e < m. This can easily understood by noticing that the potential 
term in the Dirac Hamiltonian (32) becomes negligible compared to the kinetic energy for Q — > oo. In fact, 
the introduction of the rescaled coordinate = Qxi leaves the only Q-dependence coming from the gradient 
term of the Hamiltonian. This corresponds to infinitely densely packed atoms in the model of Ref. [33] and to 
vanishing average electric potential. 

2. 1/Q — > oo: The extrema of the potential is now well-separated and one expects localized states at the minima 
(maxima) corresponding to the upper (lower) bands. Furthermore each band should be reduced to a single, 
highly degenerate energy level which corresponds to the localized states at the various minima (maxima) of 
the external potential. The semiclassical tunneling probability from a minimum to the neighbouring one is 
suppressed exponentially with increasing 1/Q, w ~ exp{—16y^mea/Q} and the problem reduces to that of the 
relativistic harmonic oscillator as far as the lowest (highest) lying states of the upper (lower) band are concerned 
[29]. 

The dispersion relations in the allowed bands alternate between convex and concave ones from band to band. If 
the Fermi-level lies inside of an allowed band, the occupied states build either a Fermi sphere or a Fermi hole in 

momentum space. For one spatial dimension the Fermi sphere is distorted to a Fermi section pi S [—kFjkp], the 
Fermi hole appears as the unoccupied interval pi G [—kF,kF] between the occupied ones, pi G [— (5/2,— /cf] and 
Pi e [kF,Q/2]. 



APPENDIX C: NUMERICAL CALCULATION OF THE ENERGY AND CHARGE DENSITY 



This appendix contains the details of the calculations of the energy and charge density. 



1. Bare expression for the energy density 



The effective Lagrangian corresponding to the action (36), 



C = -\Fp,FP^ + ^V^iT^Sp^ - \dp4,iYi> + ei^YApi, - m^V + >\{F ' F)paF'"' (CI) 



yields the energy- momentum tensor, 

dC dC - dC - 

which should be symmetrized by adding the divergence —dpf^^^^ of the third rank tensor 

= ^w^^ + L^(^jf^j'^Y - l''YY)'4}, (C3) 
8 

determined by the spin density, 
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It is easy to see that the symmetrized energy-momentum tensor is gauge invariant. The energy density operator is 



(C5) 



where the first term on the r.h.s. is the energy density of the photons and the second one is the energy density of 
the Dirac sea minus e/i multiphed by the fermion density. The last two terms appear due to the exclusion of the 
fluctuations of the collective mode. One can write £ = f (o) + £(i) + £(2) with 



1 



1 



£(1) [A, a, ^,i;,X] = ±J^ a,, {-g^'d'^F^^ + g'^°dlF°^ + ^D,!^ - Xd'J^'^ - eV^T^v) , 



(C6) 



where the lower indices indicate the powers of a. According to Eq. 40, the expectation value of the operator £ in the 
vacuum is given as £[A] = J2i=o^{i)l^^- ^® 



A/-y ^1^^ /i^y atv fitz '^e<5 c -z si?z ^olj=c=c=o 






(C7) 



with 



Zo=e^{--j-D-j-\C-G{A).Q 



(C8) 



and H£i[A) denoting the Dirac Hamiltonian (32) with the field variables replaced by functional derivatives with 
respect to the corresponding external sources, and the vacuum-to- vacuum amplitude up to the order O(e^) is 



7=1_£!/" J_ M J_ J_ J_ - J_ J_7 I 



= 1 - 



^^^^^^^ ^^^^^^^/WMVW^^^^ 




(C9) 



The 0(a°) energy, Eq. (C7), includes the energy density of the background and those of the modulated, interacting 
Dirac-sea. The 0{a) energy term is the interaction energy of the current with the fluctuations of the photon fleld, 

c rli__^/" J_ A i ^ A ^71 

" Z A,, iSr^l \6j- iSfj- iSr^y ^'^^ \5jl iSfjl ^ol^-=''=^=o 



1 

Z 




(CIO) 



The 0(0^) energy expression depends on the fluctuating field a only. Since the photon propagator in the presence of 
the background field approaches the free propagator as L ^ 00 the contribution £^2) cancels when the difference of 
the energy densities with and without the background field is considered. The vacuum to vacuum amplitude removes 
the disconnected components as expected and one finds 
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(CU) 



The propagators are calculated by means of the Lehmann representation 

na0-\^ f ^^O„-iko(xo-yo)\ f^'''^^^^0"'''^y^ , ffa'^' (a^)^^'"' (y) ] r)l^>^ - -n ^ Jk(x-y) 

^-y-2^J 2w I ko+is ^ ko-is J' 9^..2^J 2nk^+ie 

kiSi ki 

(C12) 

where f'^'^^^ (x) and gr'^i'^i (x) denote the positive and negative energy eigensolutions of the Dirac-equation (see App. A). 
The periodic background potential breaks the translational symmetry which manifests itself in changing momentum 
conservation to quasi-momentum conservation in each vertex. 



2. UV and IR divergences 

The first diagram on the r.h.s. of Eq. (Cll) represents the energy of the Dirac-sea in the presence of the background 
field, 



fsea(a, Q, ^J) = 




(C13) 



and is quadratically divergent in the absence of the background field, 



4ea(0,0,0) = - 




(C14) 



The finite, physical part of Sgea will be defined by 

41,7') (a, Q, ll) = £:,ea(a, Q, M) - 4ea(0, 0, 0), 41-') (m) = ^sea(0, 0, - .?sea(0, 0, 0) (C15) 

for the periodic and the homogeneous, normal phases, respectively. 

The convergence of £'pcr was checked numerically in the following manner. The one-loop contributions are obtained 
by taking £ of Eq. (C13) for the background field and subtracting from it the same diagram without background 
field. Let us consider first this difii'erence for vanishing chemical potential = 0, 

4i,7'^(«,g,o) = --^ ^[e(7i^(«,g,o)-47)^(o,o,o)]. (ci6) 

ki si 

By the one-by-one identification of the corresponding levels we found numerically that the magnitude jefej si {a, Q, 0) — 
Cfei si (0, 0, 0)1 is suppressed for increasing ki according to the power law 

(«, 0) - ei-X (0, 0, 0) I « kX+'^ (C17) 

with (5 « 2 > 0, cf. Fig. 10. This renders the sum absolutely convergent. The shift in the spectrum caused by the 
non-vanishing chemical potential docs not alter the UV behaviour of the sum even in the thermodynamic limit. The 
convergence of SnifJ-) = £{0, 0, fx) — £{0, 0, 0) has been checked similarly. 

IR divergences can also appear at the tadpoles where the photon line carries vanishing momentum g'^ = 0. Since 
there is actually no dynamical photon-field variable with vanishing energy and momenta, such tadpoles pose no 
problem in the homogeneous, normal vacuum [36-38]. In the periodic vacuum the photon can borrow the momentum 
nQ from the vacuum by the summation for n ^ and the tadpoles are finite. 
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3. Charge density 



In order to understand the structure of the vacuum, we need another important observable, the average charge 
density p given as the two-loop order expectation value of the operator 

/3=;^^V^7V (C18) 

The expectation value of p[A\ is taken by Eq. (40), and truncated at the two-loop order is given as 




(C19) 



These diagrams are similar to the first, fourth and fifth ones of Eq. (Cll) except that the Hamiltonian-insertion is 
replaced by a 7*^-inscrtion. The renormalization prescription (44) removes the UV divergence of the first diagram in 
Eq. (C19), too. The second diagram gives vanishing contribution for vanishing periodic background electric field due 
to Furry's theorem. The calculation of these diagrams proceeds like those for the Casimir-energy density. 



4. Two-loop diagrams 



We present now the explicit expressions for the two-loop diagrams on the r.h.s. of Eq. (41). The eigenspinors u 
and V are defined in Eq. (A3) and the diagrams containing tadpoles are 



E E ti:!^i:pv'^n:^iii + (^i - ^2)q) e "'^v^^tT^i-n + im - n2)Q), (C20) 

91 feijSi fcijSi 



e 

LT 



91 ' ' fcl, 



Vfel,Pl 
Sl,S2 



.\2 ^eni IfeS^Snz °qi-Q{ni-n2) 

fcl.Sl 
711,712 

;^klSl^^t IjPl^^A , . T ,1P1S2^0 ,.felSl,(+) X 

Z^7tl,7l2 "k"1 'K;A'^A7t2 P"2+fe"l-9l Z^71i,7t2 «"1 'afi"' 0n2 ^klSl"Pni+kn 

tpiS2 ^ '^klSl 



E 

fei,pi 

Sl,S2 



r.kisi ^.Pi^^X ,, , V ,-7PiS2^0 fcisi (-) X 

l^ni,n2^nni /kA "A7i2 "P'»2 +^"1 +9i ^"i ,>t2 "«»i Ial3^l3n2 ^kisi"Pni+k„ 

tpiS2 + f^klSl 



The exchange diagrams are given as 




= '^I'.pl'lxj dx J dyiG^^iG^^iT^,^: 



2LT 



E E 

feipigi 7li7l2 . 
S1S2 "3"4 



'^07^4 'aff"/37ii "fci,'>i+Pi."4+9i "K>^2 /kA^Atis "fei,n2-l-pi,n3-|-9i 



lai|(lai| + 4tii+41) 



(C21) 
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and 
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(C22) 



(C23) 



5. Numerical procedure 



The one-particle energy levels and spinors needed for the calculation were determined by solving numerically the 
system of linear equations (A4) and the sum over the components of the Bloch-waves was truncated for \n\ < 25. This 
procedure provided us 50 onc-particlc energy levels and spinors for each momentum in the first Brillouin zone. It was 
tested on the fermion spectrum without background field that such a truncation starts to cause noticeable error on the 
spectrum for the band index s w n when the bands are numerated in energetically increasing order. Therefore, bands 
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with s < 20 have been taken into account in the calculations of the two-loop diagrams. Such a truncation allowed 
us to detect the effects of the background field with sufficient accuracy because it was found that for ea < m the 
distortion of the dispersion relation due to the background field is only significant for states belonging to the bands 
in the vicinity of m. For about 5 bands away from m the deviation of the energy levels with and without the periodic 
background field turns practically to zero. The one- and two-loop diagrams of Eq. (41) were computed in the first 
Brillouin zone ki € [—Q/2, Q/2) at 40 and 10 points, respectively. The calculation of the one-loop diagram required 
higher numerical accuracy due to the numerical elimination of the UV divergence. We also made a test calculation 
for 20 division points which corresponded to a larger volume L and found that the numerical accuracy is about 10% 
for the two- loop contribution to the Casimir-energy density in the whole range of the parameter values. 

The amplitude a was chosen through several orders of magnitude from ea = corresponding to the perturbative 
regime to ae w m for which pair-production might occur. The wavenumber of the background field was restricted 
to be Q = 0.8, 1, 1.5 in the computation. The significantly smaller values are uninteresting from the point of view 
of the periodic state since the electrons become well localized in the limit l/Q ^ oo and it would cost too much 
energy to delocalize them. The other limit 1/Q — * is computationally time-consuming since one has to take more 
Brillouin-zones into account. Therefore the summations over the band index s, as well as over the one-particle state's 
index n have to be truncated at increasingly higher values in the formulae of App. C 1. Consequences of this very 
restricted search of the minimum of the Casimir-energy density not allowing for the higher values Q > 1.5 make 
our numerical results unreliable for large densities. The calculations were performed on an AlphaServer DS20 500 
MHz with 2 CPU-s. In particular, it took about 2 CPU hours on a single processor to compute the diagrams in Eq. 
(41) for a given set of parameters {ea,e^,/m,m) for Q = 0.8. For any given set of {Q,ea,m) we chose 50 values for 
0.5 < e/x/m < 2 in such a manner that the points between 0.98 and 1.35 were separated by 0.01. For to = 1 we took 
a = 0.0001,0.001,0.01,0.1,0.2,0.5,1.0 and Q = 0.8,1.0,1.5. For m = 0.2,0.5,2.0 and 5.0 we took only the 4 larger 
values of a. 
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FIG. 1. The scalar field 4 
compared to bs = 6e^/(/t^7r) 
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for b = 0.025, m = e = 1. 
-0.001074 from Eq. (19). 
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The slope bs = —0.001078 fitted in the central region should be 
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FIG. 3. The periodic part <pz of the scalar ground-state field configuration for b = 0.5, m = e = 1. 
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FIG. 4. Amplitude A of the periodic part (pz vs. the parameter b for various values of the electron mass m and e = 1. 
The solid line represents our perturbative estimate for the amplitude given by Eq. (24). The numerical values justify the 
perturbative estimate for small values of A. 
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FIG. 7. Density dependences of (a) the energy density, (b) the amplitude a of the photon condensate and (c) the product of 
the fiUing factor / and the wavenumber Q of the photon condensate for m = 2. The dotted hues indicate the 'critical charge 
density' pc, solid and dashed lines correspond to the periodic and normal phases, respectively. 
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FIG. 8. Dispersion relations in the first Brillouin-zone with (a) one occupied band (left) and (b) two occupied bands (right) 

sunk into the Dirac sea for ea = 0.5 and ea = 0.2, respectively, m = 2. 




FIG. 9. Allowed bands (shaded regions) for (a) undercritical ea = 0.9 (left) and (b) overcritical ea = 2 (right) static periodic 
external fields with wavenumber Q for m = e = 1. 
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FIG. 10. Momentum-dependence of the difference of the single particle energies in the expression (C16) of the one- loop 
energy for m = 2, Q = 1.1 and a = 2. 
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